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RESONANT CYCLOTRON SCATTERING IN MAGNETARS' EMISSION 
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ABSTRACT 

We present a systematic fit of a model of resonant cyclotron scattering (RCS) to the X-ray data of 
ten magnetars, including canonical and transient anomalous X-ray pulsars (AXPs), and soft gamma 
repeaters (SGRs). In this scenario, non-thermal magnetar spectra in the soft X-rays (i.e. below ~ 10 
keV) result from resonant cyclotron scattering of the thermal surface emission by hot magnetosphcric 
plasma. We find that this model can successfully account for the soft X-ray emission of magnetars, while 
using the same number of free parameters than the commonly used empirical blackbody plus power-law 
model. However, while the RCS model can alone reproduce the soft X-ray spectra of AXPs, the much 
harder spectra of SGRs below lOkeV, requires the addition of a power-law component (the latter being 
the same component responsible for their hard X-ray emission) . Although this model in its present form 
does not explain the hard X-ray emission of a few of these sources, we took this further component 
into account in our modeling not to overlook their contribution in the ~4-10keV band. We find that 
the entire class of sources is characterized by magnetospheric plasma with a density which, at resonant 
radius, is about 3 orders of magnitudes higher than uqj, the Goldreich-Julian electron density. The 
inferred values of the intervening hydrogen column densities, are also in better agreement with more 
recent estimates inferred from the fit of single X-ray edges. For the entire sample of observations, we 
find indications for a correlation between the scattering depth and the electron thermal velocity, and the 
field strength. Moreover, in most transient anomalous X-ray pulsars the outburst state is characterized 
by a relatively high surface temperature which cools down during the decay, while the properties of the 
magnetospheric electrons vary in a different way from source to source. Although the treatment of the 
magnetospheric scattering used here is only approximated, its successful application to all magnetars 
we considered shows that the RCS model is capable to catch the main features of the spectra observed 
below ~ 10 keV. 

Subject headings: radiation mechanisms: non-thermal — stars: magnetic fields — stars: neutron — 
X-rays: individual (4U 0142+614 1RXS J1708-4009 , IE 1841-045 , IE 2259+586 , 
IE 1048-5937, XTE J1810-197, IE 1547.0-5408, CXOU J1647-4552, SGR 1806-20, 
SGR 1900+14) 



1. INTRODUCTION 

The neutron star world, as we knew it until not long 
ago, appeared mainly populated by radio pulsars (PSRs, 
about 2000 objects). In the last two decades diverse, puz- 
zling classes of isolated neutron stars (NSs), with prop- 
erties much at variance with those of canonical PSRs, 
were discovered: the anomalous X-ray pulsars (AXPs) , the 
soft gamma repeaters (SGRs; Woods & Thompson 2006; 
Mereghetti 2008), the rotating radio transients (RRATs; 
McLaughlin et al. 2006), and the X-ray dim isolated neu- 
tron stars (XDINSs; Haberl 2007). Among these, the 
AXPs and SGRs are, in some sense, the most peculiar, 
since they are believed to host ultra-magnetized NSs, with 
a magnetic field « 10 14 -10 15 G, in excess of the critical 
magnetic field, B cnt = m\c? j{eh) = 4.414 x 10 13 G, at 
which the cyclotron energy equals the rest mass energy 
for an electron (Duncan & Thompson 1992; Thompson & 
Duncan 1993, 1995, 1996). 

The magnetar candidates (about fifteen known objects) 
are characterized by slow X-ray pulsations (P ~ 2-12 s) 
and large spin-down rates (P ~ 10 -10 -10 -12 fi). A dis- 



tinctive property is their high persistent X-ray luminos- 
ity (L ss 10 34 -10 36 ergs _1 ), which exceeds the spin-down 
luminosity typically, by two orders of magnitude. Thus, 
magnetar X-ray emission can not be explained in terms 
of rotational energy losses. Measurements of spin peri- 
ods and period derivatives, assuming that the latter are 
due to electromagnetic dipolar losses, lend further sup- 
port to the idea that these objects contain neutron stars 
endowed with an ultra-strong magnetic field. Although 
the magnetar model has become increasingly popular, al- 
ternative scenarios to explain the enigmatic properties of 
these sources have been proposed. Among these, models 
involving accretion from a fossil disk, formed in the super- 
nova event which gave birth to the neutron star, are still 
largely plausible (e.g. van Paradijs et al. 1995; Chatter- 
jee, Hernquist & Narayan 2000; Perna, Heyl & Hernquist 
2000). 

Magnetar X-ray emission may be qualitatively separated 
into two components, a low-energy, < 10 keV, and a high- 
energy one, > 20 keV. It is likely, although not proved 
yet, that different emission mechanisms are responsible 
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for the two components. The low energy component is 
typically fit with either a blackbody with a temperature 
kT ~ 0.3 — 0.6 keV and a power-law with a relatively 
steep photon index, T ~ 2-4, or two blackbodies with 
kTi ~ 0.3 keV and kT 2 ~ 0.7 keV (see Woods & Thomp- 
son 2006 and Mereghetti 2008 for a review). In a few 
cases the low-energy component of SCR spectra has been 
fit with a single power-law, but recent longer observations 
have shown that, also for these sources, a blackbody com- 
ponent is required (Mereghetti et al. 2005a). The high- 
energy component, discovered from four AXPs (Kuiper et 
al. 2004, 2006) and two SGRs (Mereghetti et al. 2005b; 
Molkov et al. 2005; Gotz et al. 2006) has in general a quite 
hard spectrum (modeled by a power-law), and accounts for 
about half of the bolometric luminosity of these sources. 
This makes it crucial to consider in any spectral modeling 
the whole 1-200 keV spectrum, where > 90% of the mag- 
netar emission is concentrated, instead of focussing on the 
soft X-ray range alone. Furthermore, the discovery of mag- 
netar counterparts in the radio and infrared/optical bands 
(Camilo et al. 2006; Hulleman et al. 2000) enforced the 
idea that their multi-wavelength spectral energy distribu- 
tion is by far more complex than the simple superposition 
of blackbody (BB) and power-law (PL) distributions. 

The purpose of this paper is to provide a physical inter- 
pretation of the soft X-ray component (< 10 keV) through 
a detailed analysis of magnetar spectra. Our starting point 
is the work by Thompson, Lyutikov & Kulkarni (2002, 
TLK in the following), who pointed out that resonant 
scattering in magnetar magnetospheres may explain the 
non-thermal emission observed in magnetar candidates. 
Due to the presence of hot plasma in the neutron star 
coronae, the thermal emission from the neutron star sur- 
face/atmosphere gets distorted through efficient resonant 
cyclotron scattering. Resonant cyclotron scattering has 
been first studied in the accretion columns of neutron star 
X-ray binary systems or in their atmospheres (Wasserman 
& Salpenter 1980; Nagel 1981; Lamb, Wang & Wasserman 
1990). Lyutikov & Gavriil (2006) computed, in an approx- 
imated and semi-analytical way, the effect of multiple reso- 
nant scatterings of soft photons in the magnetosphere, and 
found that the emerging spectrum is non-thermal, with 
a shape that may resemble the observed blackbody plus 
power-law. This model was preliminarily fit to the spec- 
trum of the AXP IE 1048-5937 (Lyutikov & Gavriil 2006), 
although the magnetospheric parameters were held fixed 
during the modeling. Rea et al. (2007a, b) implemented 
in XSPEC a more refined version in which also these pa- 
rameters are minimized during the fit (see ^2.2[) . and suc- 
cessfully modeled a simultaneous Swift and INTEGRAL 
observation of 4U 0142+614 . In the following, we refer to 
this XSPEC model as the RCS model, where RCS stands 
for Resonant Cyclotron Scattering. Giiver et al. (2007a, b) 
fit a similar model to two AXPs, taking into account for 
the fact that the thermal emission from the star surface 
is not a blackbody if the presence of an atmosphere is ac- 
counted for (see also More detailed, fully 3D Monte 
Carlo simulations of multiple resonant scattering in the 
star magnetosphere have been very recently presented by 
Fernandez & Thompson (2007; see also Nobili, Turolla & 
Zane 2008) but not directly applied to the data yet (this 
will be done in a subsequent paper). 



In this paper we present a systematic application of 
the RCS model to observations of all AXPs and SGRs. 
We consider the deepest X-ray pointings available up to 
now for these sources, obtained making use of the large 
throughput of the XMM-Newton satellite. For a subset 
of sources, which have been detected in the hard X-ray 
range, we also consider a joint fit with the INTEGRAL 
spectra in order to study systematically the relation be- 
tween hard and soft X-rays production mechanisms. 

The paper is organized as follows. The basic concepts 
behind the RCS model and its XSPEC implementation are 
summarized in § [5] In Sj3] we report the observations and 
the data analysis. Results of the spectral modeling are 
presented in iJH and discussed in Sj5] Conclusions follow. 

2. RESONANT CYCLOTRON SCATTERING 
2.1. The model 

Before discussing our XSPEC model and the implications 
of our results, we briefly touch on some properties of the 
RCS model which directly bear to the physical interpreta- 
tion of the fitting parameters and their comparison with 
similar parameters introduced in other theoretical models. 
The basic idea follows the original suggestion by TLK, who 
pointed out that a scattering plasma may be supplied to 
the magnetosphere by plastic deformations of the crust, 
which twist the external magnetic field and push elec- 
tric currents into the magnetosphere. The particle density 
of charge carries required to support these currents may 
largely exceed the Goldreich- Julian charge density (Gol- 
dreich & Julian 1969). Furthermore, it is expected that 
instabilities heat the plasma. 

Following this idea, Lyutikov & Gavriil (2006) studied 
how magnetospheric plasma might distort the thermal X- 
ray emission emerging from the star surface through effi- 
cient resonant cyclotron scattering. If a large volume of 
the neutron star magnetosphere is filled by a hot plasma, 
the thermal (or quasi-thermal) cooling radiation emerging 
from the star surface will experience repeated scatterings 
at the cyclotron resonance. The efficiency of the process 
is quantified by the scattering optical depth, T res , 

T res = J a res n e dl = t (1 + cos 2 a) (1) 

where 

ctt (1 + cos 2 a)uj 2 

ares = T(u-uj B )2 + ry4 (2) 

is the (non-relativistic) cross-section for electron scatter- 
ing in the magnetized regime, n e is the electrons number 
density, a is the angle between the photon propagation 
direction and the local magnetic field, T = 4e 2 ui B /3m e c 3 
is the natural width of the first cyclotron harmonic, <7t is 
the Thomson scattering cross-section, and 



3m e cujB 

Here r is the radial distance from the center of the star, 
ujb = eB/m e c is the electron cyclotron frequency, and B 
is the local value of the magnetic field. At energies cor- 
responding to soft X-ray photons, the resonant scattering 
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optical depth greatly exceeds that for Thomson scattering, 
t t ~ n e er T r, 
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This implies that even a relatively small amount of plasma 
present in the magnetosphere of the NS may considerably 
modify the emergent spectrum. 

The RCS model developed by Lyutikov & Gavriil 
(2006), and used in this investigation, is based on a sim- 
plified, ID semi-analytical treatment of resonant cyclotron 
up-scattering of soft thermal photons, under the assump- 
tion that scattering occurs in a static, non-relativistic, 
warm medium and neglecting electron recoil. The latter 
condition requires ftw <C m e c 2 . Emission from the neu- 
tron star surface is treated assuming a blackbody spec- 
trum, and that seed photons propagate in the radial di- 
rection. Magnetospheric charges are taken to have a top- 
hat velocity distribution centered at zero and extending 
up to ±/3t- Such a velocity distribution mimics a scenario 
in which the electron motion is thermal (in ID because 
charges stick to the field lines). In this respect, (3t is asso- 
ciated to the mean particle energy and hence to the tem- 
perature of the ID electron plasma. Since scatterings with 
the magnetospheric electrons occur in a thin shell of width 
H ~ Ptt/3 <C r around the "scattering sphere", one can 
treat the scattering region as a plane-parallel slab. Radi- 
ation transport is tackled by assuming that photons can 
only propagate along the slab normal, i.e. either towards 
or away from the star. Therefore, cos a = ±1 in eq. (fTJ) 
and it is r res = 2tq; the electron density is assumed to 
be constant through the slab. We notice that the model 
does not account for the bulk motion of the charges. This 
is expected since the starting point is not a self-consistent 
calculation of the currents but a prescription for the charge 
density. As a consequence, the electron velocity and the 
optical depth are independent parameters, although in a 
more detailed treatment this might not be the case (Be- 
loborodov & Thompson 2007). 

Although Thomson scattering conserves the photon en- 
ergy in the electron rest frame, the (thermal) motion of the 
charges induces a frequency shift in the observer frame. 
However, since our electron velocity distribution averages 
to zero, a photon has the same probability to undergo up 
or down-scattering. Still, a net up-scattering (and in turn 
the formation of a hard tail in the spectrum) is expected if 
the magnetic field is inhomogeneous. For a photon propa- 
gating from high to low magnetic fields, multiple resonant 
cyclotron scattering will, on average, up-scatter in energy 
the transmitted radiation, while the dispersion in energy 
decreases with optical depth (Lyutikov & Gavriil 2006). 
Photon boosting by particle thermal motion in Thomson 
limit occurs due to the spatial variation of the magnetic 
field and differs qualitatively from the (more familiar) non- 
resonant Comptonization (Kompaneets 1956). As a result, 
the emerging spectrum is non-thermal and under certain 
circumstances can be modeled with two-component spec- 
tral models consisting of a blackbody plus a power-law 
(Lyutikov & Gavriil 2006). 



2.2. The XSPEC implementation of the RCS model 

In order to implement the RCS model in XSPEC, we cre- 
ated a grid of spectral models for a set of values of the 
three parameters (3t, T res and T. The parameter ranges 
are 0.1 < (3t < 0.5 (step 0.1; @t is the thermal veloc- 
ity in units of c), 1 < r res < 10 (step 1; tau res is the 
optical depth) and 0.1 keV< T < 1.3 keV (step 0.2 keV; 
T is the temperature of the seed thermal surface radia- 
tion, assumed to be a blackbody). For each model, the 
spectrum was computed in the energy range 0.01-10 keV 
(bin width 0.05 keV). The final XSPEC atable spectral 
model has therefore three parameters, plus the normal- 
ization constant, which are simultaneously varied during 
the spectral fitting following the standard x 2 minimiza- 
tion technique. In Fig. we show the comparison between 
a blackbody model and our RCS model. We stress again 
that our model has the same number of free parameters 
(three plus the normalization) than the blackbody plus 
power-law or two blackbody models (/3t, i~ res , T, plus the 
normalization, compared to kT, T (or kT2), plus two nor- 
malizations); it has then the same statistical significance. 
We perform in the following section a quantitative compar- 
ison between the RCS model and other models commonly 
used in the soft X-ray range. However, note that here the 
RCS model is meant to model spectra in the 0.1-10 keV 
energy range. For all sources with strong emission above 
~ 20 keV, the spectrum was modeled by adding to the 
RCS a power-law meant to reproduce the hard tail (see 
§2] for details). This power-law does not have (yet) a clear 
physical meaning in our treatment, but since it contributes 
also to the 0.1-10 keV band, our RCS parameters depend 
on the correct inclusion of this further component. 

3. OBSERVATIONS AND DATA ANALYSIS 

Before discussing our data analysis, we would like to 
outline the choices we made in selecting the datasets to 
be used in this work. Aim of this paper is to show how 
the RCS model can account for the X-ray spectra of both 
steady and variable AXPs and SGRs. Detailed spectral 
modeling requires high-quality data and this led us to con- 
sider only the highest signal-to-noise ratio datasets avail- 
able to date for these sources. We selected then only 
those magnetar candidates having XMM-Newton spec- 
tra with a number of counts > 10 5 and did not include 
short (< 10 ks) XMM-Newton exposure^, Chandra or 
Swift observations. Fortunately most of the magnetars 
met the above criterion, but our choice resulted in the ex- 
clusion of CXOUJ0100-7211, AXJ1844-0258, SGR 0526- 
66, SGR 1627-41 and SGR 1801-23; they are no further 
considered in the present investigation^. The remaining 
sources are divided into three groups, as follows. 

• A set of AXPs which emit in the hard X-ray range, 
and also happen to be "steady" emitters or show- 
ing moderate flux and spectral variability (flux 
changes less than a factor of 5; with the exception 
of IE 2259+586 , see also below). These long-term 
changes are not considered in the following (see ij l4.ll 



for details). This group comprises: 4U 0142+614, 

6 Except for IE 1841-045 , for which only a single short XMM-Newton observation is available. 

7 While this paper was approaching completion, Tiengo, Esposito & Mereghetti (2008, ApJ submitted) reported a detailed 0.1-10 keV spectrum 
for CXOU J0100-7211 . In their paper, the successful application of our RCS model to this source is presented. 
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1RXS J1708-4009, IE 1841-045, and IE 2259+586 . 
When more than one XMM-Newton observation 
was available, we chose the dataset with the longest 
exposure time and least affected by background 
flares. 

• A set of "transient" AXPs (often labeled TAXPs), 
which includes XTE J1810-197, IE 1547.0-5408, 
and CXOUJ1647-4552. To these we add IE 1048- 
5937 , in the light of the recent detection of large 
outbursts from this source (Mereghetti et al. 2004; 
Gavriil et al. 2006; Tarn et al. 2007; Campana 
& Israel 2007), and of its spectral similarities with 
canonical TAXPs. In order to follow the spectral 
evolution without being encumbered with unneces- 
sary details, we selected only three XMM-Newton 
spectra for each source, also when more observa- 
tions were available (e.g. for IE 1048-5937 and 
XTE J1810-197). The three chosen datasets cor- 
respond to the two most diverse spectra and to an 
"intermediate" state. 

• A set of SGRs, which comprises SGR 1806- 
20 (three observations covering epochs before and 
after the giant flare of 2004 December 27), and 
SGR 1900+14. 

For all the sources in the first group (except 
IE 2259+586) and for SGR 1900+14 we also considered 
INTEGRAL data. Although INTEGRAL and XMM- 
Newton observations were not always simultaneous, the 
absence of large spectral variability in these sources jus- 
tifies our choice. In particular, for SGR 1900+14 care 
has been taken to select data within periods in which the 
source was relatively steady. Although AXP IE 2259+586 
and SGR 1806-20 have been also detected above 20keV 
(Kuiper et al. 2006; Mereghetti et al. 2005b; Molkov et 
al. 2005), the INTEGRAL X-ray counterpart of the for- 
mer is too faint to extract a reliable spectrum, while the 
highly variable hard and soft X-ray spectrum of the latter, 
together with the non simultaneity of the XMM-Newton 
and INTEGRAL observations, would make any attempt 
to model its 1-200 keV spectral energy distribution mean- 
ingless. 

The following subsections provide some details on the 
observations and data analysis; a comprehensive log, with 
the exposure times and epochs of each observation, is pro- 
vided in Tab.[U 

3.1. XMM-Newton: soft X-rays 

All soft X-ray spectra were collected by the XMM- 
Newton EPIC-pn instrument (Jansen et al. 2001; Striider 
et al. 2001), which has the largest sensitivity in the 1- 
10 keV band. In order to have a homogeneous sample 
of spectra, we re-analysed all the data using the latest 
SAS release 7.1.0. We employed the most updated cali- 
bration files available at the time the reduction was per- 
formed (August 2007). Standard data screening criteria 
(e.g. cleaning for background flares) were applied in the 
extraction of scientific products. We used FLAG= and 
PATTERN between 0-4 (i.e. single and double events) 
for all the spectra. We have checked that spectra gener- 
ated with only single events (i.e. PATTERN^ 0) agreed 



(apart from normalization factors) with those generated 
from single and double events. All the EPIC-pn spectra 
were rebinned before fitting, using at least 30 counts per 
bin and not oversampling the resolution by more than a 
factor of 3 (see Rea et al. 2005, 2007c for further details 
on our XMM-Newton data analysis and reduction). 

3.2. INTEGRAL: hard X-rays 

In order to take into account in our spectral 
modeling the contribution of the hard X-ray emis- 
sion of 4U 0142+614, 1RXSJ1708-4009, IE 1841-045 and 
SGR 1900+14, we used the hard X-ray spectra derived 
from INTEGRAL data. We selected and analyzed all 
publicly available IBIS (Ubertini et al. 2003) pointings, 
making use of ISGRI (Lebrun et al. 2003), the IBIS low 
energy detector array working in the 15keV— lMeV en- 
ergy range. Data were collected for all pointings within 
12° from the direction of each source, for a total 2544, 
1351, 1894 and 1535 pointings of 2-3 ks each, for the three 
AXPs and the SGR, respectively. Given the low hard X- 
ray flux of these sources, we added all the data in order to 
have statistically significant detections. 

We processed the data using the Offline Scientific Anal- 
ysis (OSA) software provided by the INTEGRAL Science 
Data Centre (ISDC) v6.0. We produced the sky images of 
each pointing in 10 energy bands between 20 and 300 keV, 
and added them in order to produce a mosaicked image. 
Due to the faintness of the sources we could not derive 
their spectra from the individual pointings, so following 
e.g. Gotz et al. (2007), we used the count rates of the 
mosaicked images to build the time averaged spectrum of 
each source. 

4. SPECTRAL ANALYSIS AND RESULTS 

All the fits have been performed using XSPEC version 
11.3 and 12.0, for a consistency check. A 2% systematic 
error was added to the data to partially account for uncer- 
tainties in instrumental calibrations. A constant function 
has been fitted when using both XMM-Newton and IN- 
TEGRAL data to account for inter-calibration uncertain- 
ties (the values of the constant in the Tables are relative to 
XMM-Newton set to unity). The 0.5-1 keV energy range 
was excluded from our spectral fitting because: i) this is 
the band where most of the calibration issues lay (Haberl 
et al. 2004), and ii) emission in this energy range is mostly 
affected by interstellar absorption, and by the choice of the 
assumed solar abundances. Given the high column density 
of all magnetars, and the large uncertainties in the abun- 
dances (probably not even solar) in their directions, this 
may lead to spurious features. We checked that for all 
our targets, the values of Nh derived fitting the 1-10 keV 
EPIC-pn spectra, are consistent (within the errors) with 
those obtained using the 0.5-10 keV range in the same data 
set. We notice that the absorption value derived here for 
the blackbody plus power-law or two blackbodies models 
is, on average, slightly higher than that reported in the lit- 
erature for the same model. This is due to our choice of us- 
ing the more updated solar abundances by Lodders (2003), 
instead of the older ones from Anders & Grevesse (1989). 
This does not affect the other spectral parameters, which 
are in fact consistent with those previously published for 
the same data sets. For all the fits we used photoelectric 
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cross-sections derived from Balucinska-Church & McCam- 
mon (1992). 

We raise the caveat that no attempt has been made here 
to distinguish the pulsed from the non-pulsed emission of 
these objects, and to model the spectral variability with 
phase observed in most of these sources. This will be the 
subject of a future investigation. 

4.1. AXPs: the hard X-ray emitters 

In this section we first consider the AXPs with de- 
tected hard X-ray emission, which also coincides with 
the marginally variable AXPs, with the exception of 
IE 2259+586 (Kaspi et al. 2003; Woods et al. 2004; see be- 
low). We recall that, strictly speaking, these hard X-ray 
emitting AXPs are not "steady" X-ray emitters. Subtle 
flux and spectral variability was discovered in 1RXS J1708- 
4009 and 4U 0142+614. In particular, 1RXS J1708-4009 
showed a long term, correlated intensity-hardness variabil- 
ity (both in the soft and hard X-rays), most probably re- 
lated to its glitching activity (Rea et al. 2005; Campana 
et al. 2007; Gotz et al. 2007; Dib et al. 2007a; Israel et 
al. 2007a). 4U 0142+614 showed a flux increase of ~ 10% 
(also correlated with a spectral hardening) following the 
discovery of its bursting activity (Dib et al. 2007b; Gonza- 
lez et al. 2007). Furthermore, thanks to a large KXTE 
monitoring campaign, long-term spin period variations 
and glitches were discovered in 4U 0142+614 1RXS J1708- 
4009 , and IE 1841-045 , i.e. the three AXPs which are the 
brightest both in the soft and hard X-ray bands (Gavriil 
& Kaspi 2002; DalPOsso et al. 2003; Dib et al. 2007a; 
Israel et al. 2007a). 

Since these flux variations are rather small, we have cho- 
sen to model only the XMM-Newton observation closest 
to the INTEGRAL one (for 1RXS J1708-4009 only one 
XMM-Newton observation is available though). Our re- 
sults from the spectral modeling of the 1-200 keV spec- 
trum of 4U 0142+614, 1RXSJ1708-4009, and IE 1841- 
045 are summarized in Table and shown in Fig. O 

The case of IE 2259+586 is rather different: it showed 
a large outburst (more than one order of magnitude flux 
increase) detected by i?XTE , during which also bursting 
activity was detected (Kaspi et al. 2003). However, in 
the XMM-Newton observations pre and post outburst, 
the source showed fluxes which differ only by a factor of 
3 (Woods et al. 2004). Furthermore, it was observed to 
emit up to ~30 keV by the HEXTE instrument on board 
of i?XTE (Kuiper et al. 2006) and by INTEGRAL , but 
unfortunately it is too faint in the latter observation to 
extract a spectrum. We then decided to model only the 
deepest XMM-Newton observation taking into account of 
the >10keV component by adding a power-law with pho- 
ton index fixed at the HEXTE value (Kuiper et al. 2006). 
This is because sizable residuals are present at the highest 
energies when the XMM-Newton spectrum is modeled ei- 
ther with the BB+PL or the RCS model. A satisfactory 
fit requires, in both cases, the addition of a hard X-ray 
power-law component (see also Tableland Fig. [2]). 

Summarizing, the only source that can be considered (so 
far) a genuine "steady" X-ray emitter, among the AXPs 
with hard X-ray emission is IE 1841-045 . It is interesting 

8 This is because the absorption model tends to increase the Nh value 
eventually diverges approaching E=0. 



to note that this is also the only AXP for which a black- 
body plus a single power-law reproduces well the entire 
1-200 keV spectrum, while for the other hard X-ray emit- 
ting AXPs two power-laws are required. In this respect, 
the spectral distribution of IE 1841-045 resembles the one 
of the SGRs (see also 

In all cases we found that Nh, as derived from the RCS 
model, is lower than (or consistent with) that inferred 
from the BB+2PL fit (or BB+PL in the case of IE 1841- 
045 ) , and consistent with what derived from fitting the 
single X-ray edges of 4U 0142+614, IE 2259+586, and 
1RXS J 1 708-4009 (Durant & van Kerkwijk 2006). This is 
not surprising, since the power-law usually fitted to mag- 
netar spectra in the soft X-ray range is well known to 
cause an overestimate in the column densitj0 The surface 
temperature we derived fitting the RCS model is system- 
atically lower than the corresponding BB temperature in 
the BB+2PL or BB+PL models, and is consistent with 
being the same (~ 0.33 keV) in the four sources. On the 
other hand the thermal electron velocity and the optical 
depth are in the ranges 0.2-0.4 and 1.0-2.1, respectively. 
Concerning the hard X-ray power-law, we find that the 
photon index is, within the errors, the same when fitting 
the RCS or the BB+2PL or BB+PL models (note that for 
IE 2259+586 it was kept fixed), while the hard PL nor- 
malization is larger in the RCS case with respect to the 
BB+2PL model. Both the soft and the hard X-ray fluxes 
of all these AXPs derived from the RCS fitting are consis- 
tent with those implied by the usual BB+2PL fitting. 

4.2. AXPs: the "transients" 

"Transient" AXPs have been discovered only very re- 
cently, when an increase in the X-ray flux by a factor 
~ 100 over the value measured a few years before was 
observed in XTEJ1810-197 (Ibrahim et al. 2004; Got- 
thelf et al. 2004). Later on, new TAXPs have been ob- 
served showing large flux and spectral variations, e.g. 
CXOU J1647-4552 (Muno et al. 2007) and IE 1547.0-5408 
(Gelfand & Gaensler 2007; Camilo et al. 2007a; Halpern et 
al. 2007). Very intriguing is the discovery of pulsed radio 
emission correlated with the outbursts of XTE J1810-197 
and IE 1547.0-5408 (Camilo et al. 2006, 2007a), while so 
far only upper limits have been set on the radio emission 
from CXOUJ1647-4552, IE 1048-5937 and other AXPs 
(Burgay et al. 2006, 2007; Camilo et al. 2007b). 

It is not clear whether AXPs and TAXPs are indeed 
two distinct groups of sources. During the past few years 
it has became increasingly evident that flux variations of 
different magnitudes also occur in "steady" AXPs, pos- 
sibly related to their bursting and glitching activity (see 
q4.1[) . Furthermore, bursts have been observed also dur- 
ing the outbursts of the TAXP XTEJ1810-197 (Woods 
et al. 2005) and CXOU J1647-4552 (Muno et al. 2007), 
the latter also showing a large glitch (Israel et al. 2007b). 
However, in this paper we maintain the distinction be- 
tween TAXPs and AXPs, partly for historical reasons, and 
partly because the two classes may indeed have different 
spectral properties, with the TAXPs being characterized 
by much softer X-ray spectra, and by the lack, so far, of 
detection at energies > lOkeV. 

in response of the steep rise of the power-law at low energies, which 
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The results of the TAXPs spectral modeling are sum- 
marized in Tables 4, 5, 6, 7, and shown in Figs. [4] El 
[5J El Also in this case, we chose to model up to three 
spectra representative of the flux and spectral variabil- 
ity of these sources. Again, Nh derived with the RCS 
model is lower than (or consistent with) that inferred 
from the more common BB+BB fitting for XTEJ1810- 
197, CXOUJ1647-4552, and IE 1547.0-5408 , and signifi- 
cantly lower in the case of the BB+PL model applied to 
IE 1048-5937 (and consistent with that derived by Du- 
rant & van Kerkwijk 2006). We also found that the RCS 
model can easily account for all the spectral and intensity 
changes in the TAXPs. With the exception of XTE J1810- 
197 , the surface temperature we derive for all the TAXPs 
is lower than, or consistent with, that of the blackbody 
in the BB+PL or BB+BB model (for the BB+BB model, 
we refer to the BB with the lowest temperature). How- 
ever, considering only the RCS model, it is evident for 
XTEJ1810-197, IE 1547.0-5408, and CXOU J1647-4552 
that the outburst state has a high surface temperature 
which cools down during the decay, while for IE 1048-5937 
this trend is less clear. Furthermore, for all the TAXP but 
CXOU J1647-4552, (3t increases during the outburst de- 
cay. The behavior of T res is less homogeneous: this param- 
eter decreases with decaying flux in XTE J1810-197and 
IE 1048-5937, remains qconstant in IE 1547.0-5408 , and 
shows an increase during the outburst decay in the case of 
CXOU J1647-4552. Also for these transient sources, the 
fluxes derived by the empirical model and the RCS model 
are consistent. 

4.3. SGRs 

Finally, we consider the 1-10 keV and 1-200 keV 
emission of SGR 1806-20 (see Table 8 and Fig. 8) and 
SCR 1900+14 (see Tableland Fig.©, respectively. It has 
been already noticed that the hard X-ray emission of SGRs 
is quite different from that of AXPs (see fc|4.1|) . In fact, 
with the exception of IE 1841-045 , the spectra of AXPs 
show a clear turnover between 10 and 20 keV (see Fig. ^ 
and the fit requires an additional spectral component. In- 
stead, the hard X-ray emission of SGRs seems the natural 
continuation of the non-thermal component which is dom- 
inant in the 1-10 keV energy range. This is why we can 
use a BB (or RCS) plus a single power-law in the entire 1- 
200 keV range for SGR 1900+14 , while for the hard X-ray 
emitting AXPs we were forced to add a second power-law 
to the BB+PL model. 

Similar considerations hold for SGR 1806-20 , in which 
case we model the 1-10 keV emission by adding a power- 
law component which is intended to account for the contri- 
bution of the hard X-ray emission in the soft X-ray range. 
For the latter SGR we modeled three X-ray observations 
taken before and after the Giant Flare of 2004 Decem- 
ber 27 (Hurley et al. 2005; Palmer et al. 2005). We 
found that the Njj value is consistent within 1 a between 
the BB+PL and the RCS+PL models, and the power-law 
contribution and the photon index vary among the three 
spectra in a similar fashion for the two models. Also, in 
the RCS+PL model the surface temperature remains con- 
stant within the errors until before the Giant Flare, and 
then becomes very low after one year. Besides the temper- 
ature, the spectral variability is accounted for by changes 



in the parameters describing the magnetospheric currents, 
with (3t and r res varying in the ranges 0.14-0.5 and in the 
2.2-4.3, respectively. 

In the SGR 1900+14 1-200 keV spectrum, we found 
consistent Nh and spectral index values between the 
BB+PL and RCS+PL models, and a RCS surface temper- 
ature significantly lower than the corresponding BB tem- 
perature. In all the SGR observations, the derived fluxes 
are consistent among the two models. 

5. DISCUSSION 

Before discussing our results and the physics we can de- 
rive from our model, we would like to stress once again that 
the RCS model involves a number of simplifications (see 
i j2.ip . One is the assumption of a single temperature sur- 
face emission. Current-carrying charges will hit and heat 
the star surface, generally inhomogencously (TLK). In ad- 
dition, the emission emerging from the surface is likely to 
be non-Plankian. While the presence of an atmosphere 
on top the crust of a magnetar remains a possibility (see 
Giiver et al. 2007a, b), its properties, are then likely differ- 
ent from those of a standard (in radiative and hydrostatic 
equilibrium) atmosphere on, e.g., a canonical isolated cool- 
ing neutron star (see e.g. Ho & Lai 2003; van Adelsberg & 
Lai 2006). The extreme field and (relatively) low surface 
temperature (< 0.5 keV) of magnetar candidates may also 
be suggestive of a condensed surface, at least if the chem- 
ical composition is mainly Fe (see Turolla, Zane & Drake 
2004). In the light of these considerations, and in the ab- 
sence of a detailed model for the surface emission, and for 
the atmosphere of strongly magnetized NSs constantly hit 
by returning currents, we restricted ourself to a blackbody 
approximation for the seed thermal photons. 

In spite of these simplifications, we find that the RCS 
model can describe the soft X-ray portion of the whole 
set of magnetar spectra we have considered, including the 
TAXPs variability, by using only three free parameters 
(plus a normalization factor). This is the same number of 
degrees of freedom required by the blackbody plus power 
law model, commonly used to fit this energy band. 

5.1. Magnetar magnetospheric properties 

One of the most interesting outcomes of our analysis is 
the measure of the magnetospheric properties of magne- 
tars. In all sources, steady and variable ones, the value of 
r res is in the range of ~ 1-6. This suggests that the en- 
tire class of sources are characterized by similar properties 
of scattering electrons, their density and their (thermal) 
velocity spread. An optical depth tq = r res /2 requires a 
particle density n e (see eq. [5] ) which can be easily inferred 
considering: 

T « 1.8 x 10- 2 Vr sc , (5) 

where r sc is the radius of the scattering sphere 

nr, ( B \ 1/3 /lkeV\ 1/3 

r sc ^8R NS -— , (6) 

\Bcrit J \ hW B J 

Rns is the neutron star radius and B crit ~ 4.4 x 10 13 G 
is the quantum critical field. By taking a typical photon 
energy of ~ 1 keV, Rns ~ 10 6 cm and B ~ l0B cr it, 
we get n e w 1.5 x 10 13 r res cm -3 . This is several or- 
ders of magnitude larger than the Goldreich-Julian den- 
sity (Goldreich & Julian 1969) at the same distance, 
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no j ~ n e irr sc /($T res Ri c ) — 2 x 10 10 cm 3 (where R tc is 
the light cylinder radius and we took P ~ 10 s). While the 
charge density is large when compared with the minimal 
Goldreich- Julian density, it provides a negligible optical 
depth to non-resonant Thomson scattering. Only the reso- 
nant cyclotron scattering makes an efficient photon boost- 
ing possible. 

Our present model does not include a proper treatment 
of magnetospheric currents, so that r res is a free parameter 
related to the electron density. Nevertheless, it is useful 
to compare the values of the optical depth inferred here 
to those expected when a current flow arises because a 
steady twist is implanted in the star magnetosphere, as 
in the case investigated by TLK under the assumption 
of axysimmetry and self-similarity If the scattering par- 
ticles have a collective motion (bulk velocity fibuik), the 
efficiency of the scattering process is related to T res f3b u ik 
(e.g. Nobili, Turolla & Zampieri 1993). This quantity is 
shown as a function of the magnetic colatitude in Fig. 5 of 
TLK for different values of the twist angle, A<f>iy_s- By 
assuming fibuik = 1 and integrating over the angle, we get 
the average value of the scattering depth as a function of 
A</>jv_s, which is shown in Fig. [TO] The curves correspond- 
ing to a different value of ftbuik can be obtained simply by 
reading the quantity shown in Fig. [10] as T res (3b u ik and by 
rescaling the j/-axis. As we can see, a value of r res ~ 1 
is only compatible with very large values of the twist an- 
gle (i.e A0jv_s > 3), while typical values of r res ~ 2, as 
those obtained from some of our fits, require Pbuik ^ 0.5 
to be compatible with A(f>N-s ~ 3 (the smaller is Pbuik, 
the smaller is the value of the twist angle). This is consis- 
tent with the fact that the RCS model has been computed 
under the assumption of vanishing bulk velocity for the 
magnetospheric currents, and it is compatible with TLK 
model only when in the latter it is ftbuik -C 1. 

5.2. Comparison between AXPs and SGRs 

In the last few years the detection of bursts from AXPs 
(Gavriil et al. 2002; Kaspi et al. 2003) strengthened their 
connection with SGRs. However, the latter behave differ- 
ently in many respects. Below ~ 4 keV, the SGRs emission 
can be described either by a blackbody or an RCS compo- 
nent. At higher energies though (> 4keV), their spectra 
require the addition of a power-law component, which well 
describes the spectrum until ~ 200 keV. The non-thermal 
component dominates their spectra to the point that the 
choice of a blackbody or the RCS model at lower energies 
does not affect significantly the value of the hard X-ray 
power-law index, nor the energy at which this component 
starts to dominate the spectrum (see e.g Tab. 9 and Fig. 9). 
The spectra of SGRs are then strongly non-thermally dom- 
inated in the 4-200 keV range. 

The case of the AXPs is different (with the excep- 
tion of IE 1841-045, see below). These sources show 
a more complex spectrum, with an evident non-thermal 
component below ~ 10 keV, apparently different from 
that observed at higher energies. For the AXPs detected 
at energies >20keV, the spectrum can be described by 
a RCS component until 5-8 keV, above which the non- 
thermal hard X-ray component becomes important, and 
(e.g. for 1RXS J1708-4009 and 4U 0142+614) dominates 
until ~ 200 keV. In the case of the F3B+2PL model in- 



stead, the non-thermal component responsible for the hard 
X-ray part of the spectrum starts to dominate only above 
~ 10 keV (see e.g. Figs. [5] and [3]). This is important, 
because the measurement of a down-break of the hard X- 
ray power-law has remarkable physical implications and 
may prove useful in constraining the physical parameters 
of the model for the hard X-ray emission. It is worth not- 
ing that the photon index of the hard X-ray component 
in AXPs does not strongly depend on the modeling of the 
spectrum below lOkeV, while, its normalization and, as 
a consequence, the value at which the hard tail starts to 
dominate the spectrum, do. 

In this picture IE 1841-045 seems an exception. From 
the spectral point of view, IE 1841-045 appears as the 
more SGR-like among the AXPs. Its multi-band spec- 
trum can be well fitted by a BB+PL or RCS+PL model, 
with parameters very similar to SGRs (compare Figs. [91 
and Tables [5]). This may suggest that this source 
is a potential transition object between the two classes. 
However, at variance with the SGRs, this source seems to 
be the least active bursters among AXPs. Note that, at 
variance with the other magnetars, in the case of the two 
SGRs and IE 1841-045 , our model requires two additional 
free parameters, with respect to the BB+PL, to account 
for the hard X-ray power-law. 

The fact that hard X-ray spectra detected from AXPs 
are much flatter than those of SGRs may also suggest a 
possible difference in the physical mechanism that powers 
the hard tail in the two classes of sources. Within the mag- 
netar scenario, Thompson & Beloborodov (2005) discussed 
how soft 7-rays may be produced in a twisted magneto- 
sphere, proposing two different pictures: either thermal 
bremsstrahlung emission from the surface region heated by 
returning currents, or synchrotron emission from pairs cre- 
ated higher up (~ 100 km) in the magnetosphere. More- 
over, a third scenario involving resonant magnetic Comp- 
ton up-scattering of soft X-ray photons by a non-thermal 
population of highly relativistic electrons has been pro- 
posed by Baring & Harding (2007). It is interesting to note 
that 3D Monte Carlo simulations (Fernandez & Thomp- 
son 2007; Nobili, Turolla & Zane 2008) show that multiple 
peaks may appear in the spectrum. In particular, in the 
model by Nobili, Turolla & Zane (2008), a second "hump" 
may be present when up-scattering is so efficient that pho- 
tons start to fill the Wien peak at the typical energy of the 
scattering electrons. The change in the spectral slope may 
be due, in this scenario, to the peculiar, "double-humped" 
shape of the continuum. The precise localization of the 
down-break is therefore of great potential importance and 
might provide useful information on the underlying phys- 
ical mechanism responsible for the hard emission. 

The RCS model applied to the evolution of the outbursts 
of the TAXPs known up to now shows how the outburst 
may results from a heating of the NS surface, which slowly 
cools in a timescale of months/years. AXPs outbursts are 
thought to be caused by large scale rearrangement of the 
surface/magnetospheric field, either accompanied or trig- 
gered by fracturing of the NS crust. It is worth noticing 
that from our modeling we find that the surface tempera- 
ture cools down during the outburst decay, while the mag- 
netospheric characteristics change in a different way from 
source to source. 
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5.3. Correlations 

The quite large number of observations we analyzed 
(both relative to different sources and to single sources 
in different emission states) allows to search for possible 
correlations among the various quantities, both in the en- 
tire sample, i.e. looking at the population of magnetar 
candidates at large, and in the time evolution of a single 
source. 

Fig. 11 summarizes the results of our spectral fits. The 
various panels show how the three model parameters (T, 
T res and Pt) are related to the X-ray luminosity in the 
1-10 keV band (£i-iokev) and to the magnetic field B. 
The latter is derived from P and P, assuming that the 
magnetic field is a core-centered dipole and the spin-down 
is due magnetic dipole radiation. 

An inspection of the panels in Fig. 11 does not reveal 
any obvious correlation for the entire set of observations. 
To verify this, we have run a Spearman rank test and we 
only found a positive correlation between B and both r res 
and Pt (deviation from the null hypothesis at about the 
93% and 89% confidence level, respectively). No correla- 
tions with a significance level above ~ 65% were found in 
all the other cases. Both parameters @t and r res control 
the scattering efficiency, but the meaning of their correla- 
tion with the field strength, which seems to be direct in the 
case of the optical depth and inverse in the case of the ther- 
mal velocity (Fig. 11) is not of immediate interpretation. 
The optical depth scales as n e r/B (see eq. [3]). If we make 
again a comparison with the twisted magnetosphere model 
(TLK), in which n e oc B/r, this is not expected. Taken 
face value, an increase of the optical depth with increasing 
B implies that the product n e r grows more rapidly than 
B. Since both in the RCS model and in TLK the scat- 
tering radius is oc B 1 ^ 3 , this implies that n e should grow 
faster than what expected in a self-similar magnetostatic 
configuration. Furthermore, it is interesting to note 

On the other, we caveat that these considerations are 
largely model dependent and, in order to assess this is- 
sue, a detailed treatment of the magnetosphere, including 
more realistic profiles for the electron density and velocity 
distribution, is needed. 

As discussed earlier, a more interesting trend is found 
restricting to observations of the same source at differ- 
ent epochs. In many transient AXPs (e.g. XTEJ1810- 
197, IE 1547.0-5408, and CXOU J1647-4552) we observe 
a clear correlation between the surface temperature and 
the X-ray luminosity, which is expected since in the RCS 
model an enhanced surface thermal emission produces 
more seeds for resonant up-scattering. However, once 
again there is no clear trend relating changes in T res and (3t 
to changes in luminosity for the entire TAXP sample. In 
most transient sources at least one of these two parameters 
increases with flux, and this may be enough to guarantee 
that the spectrum hardens at larger luminosities, but in no 
case there is a simultaneous increase or decrease of both 



r res and Pt during the outburst decay. Whether this is 
due to a degeneracy in the model parameter space or it 
reflects a real trend is not clear at present. 

6. CONCLUSION 

In this paper we showed that the soft X-ray emission 
of magnetars can be explained by resonant cyclotron scat- 
tering of their thermal surface emission by a cloud of hot 
magnetospheric electrons. This model satisfactorily repro- 
duces the spectral shape of all magnetars soft X-ray emis- 
sion, using the same number of free parameters than the 
widely used blackbody plus power-law model (except for 
the SGRs where the much harder spectrum below lOkeV, 
still requires the addition of a power-law on top of the res- 
onant cyclotron scattering model, being the same power- 
law component responsible for their hard X-ray emission). 
This means that the RCS model not only catches the main 
features of the thermal and non-thermal components ob- 
served in these sources below ~ 10 keV, but also success- 
fully provides a quantitative interpretation. For the mag- 
netars presenting an hard X-ray emission we included this 
further component in order to take into account in our 
modeling of the contribution of this component down to 
the soft X-ray part of the spectrum. 

This work represents one of the first attempts to infer 
some physical values from the 1 — 10 keV spectra of mag- 
netars. Future refinements are in progress, in order to 
improve the RCS model from a ID analytical model to- 
ward a 3D Monte Carlo based code (as the more advanced 
codes developed by Fernandez & Thompson 2007 and No- 
bili, Turolla & Zane 2008). Furthermore, this model even- 
tually applied to the detailed spectra that XEUS and/or 
Con-X will possibly make available in the near future, ap- 
pear a promising step toward the complete understanding 
of the physics behind magnetars soft X-ray emission. 
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Table 1 

Log Of The XMM-Newton and INTEGRAL Observations Analysed In This Paper. 

XMM-Newton 



Source Date (YYYY/MM/DD) Exposure (ks) 



/ITT 01 A A 


900/1 /O^ /01 
ZUU4/ Uo/ Ul 


A A 
44 


T170& zLOOQ 


900^ /OR /9R 


AJ\ 

4<J 


1 V, 1 841 04^ 


9009/1 /07 


U 


1 E 2259+586 


2002/06/1 1 


52 


IE 1048-5937 


2003/06/16 


69 




2005/06/17 


32 




2007/06/14 


48 


XTEJ1810-197 


2004/09/18 


28 




2005/09/20 


42 




2006/03/13 


51 


IE 1547.0-5408 


2006/08/21 


47 




2007/08/09 


16 


CXOUJ1647-4552 


2006/09/16 


80 




2006/09/22 


20 


SGR 1806-20 


2003/04/03 


55 




2004/10/06 


19 




2005/10/04 


33 


SGR 1900+14 


2005/09/17 


30 


INTEGRAL 


Source 


Date (YYYY/MM/DD) Exp 


3sure (Ms) 


4U0142+614 


2003/03/03-2006/08/13 


1.9 


1RXS J 1708-4009 


2003/02/28-2005/10/02 


2.7 


IE 1841-045 


2003/03/10-2006/04/28 


4.0 


SGR 1900+14 


2003/03/06-2006/09/26 


3.7 







kT = 0.8keV 




kT = 0.2keV 








\\ 
"\ \ 






\ \ 









12 5 10 



Energy (keV) 

Fig. 1. — Distorsion of a seed blackbody spectrum through resonant cyclotron scattering onto magnctosfcric electrons, for two values of the 
blackbody temperature, 0.2 keV and 0.8 keV. Black lines show the RCS model for (if = 0.2 and T res = 2, 4, 8 (from bottom to top), while 
grey lines are relative to (3t = 0.4 and r res = 2, 4, 8 (from bottom to top). The normalizations of the various curves are arbitrary. 
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Table 2 

Spectral Parameters: 4U 0142+614, 1RXS J1708-4009 , and IE 1841-045 



AXPs 



Parameters 



4U 0142+614* 



1RXSJ1708-4009* 



IE 1841-045 



BB+2PL 

-L- ' -0.02 

1.01 



RCS+PL 
n o-i+u.ub 

1.10 



BB+2PL 

I 1+0. Ob 

1 - yi -0.06 

1.05 



RCS+PL 

1 67+u.ub 
-■-•"'-o.os 

0.80 



BB+PL 

5+0.4 



RCS+PL 

9 r 7 +U.13 
z - ' -0.15 

1.09 



N ff 

constant 



2.38~ 



1.02 



o.i 



kT (keV) O.i:! :;;,!;; 0.30±g:8g 0.47lg;gi imtfcgg 0.5lig : g 0.39±g:8g 

BB norm 8.7lg"| x 10~ 4 2.4±g-^ x 10" 4 2 - 4 ±o"i x 10 ~ 4 

r 4 14+0.04 2 7n+ 08 

x l 4 - 1 *-0.04 z -' u -0.08 

PLi norm 0.30±°;°1 0.016±8;gg| 



a n oq+0.05 n qo+0.03 n oo +0.05 

PT 0.33_ 005 U.dH_ 0-03 U.26_ 005 

1 q+0.2 9 1+0.2 1 10+0.3 

'res — 0.2 *■■>-— 0.2 x - lo -0.2 

RCS norm 4.5±g;| x 10" 3 8.1±i:| X 10~ 4 3.1±?;f x 10~ 4 



r n 7S+ - 1 i 1+ - 1 n 7R+0.1 i n+o- 1 1 47+0.04 , 17+0.05 

L 2 U+8_ 07 El_ 01 U.7b_ 01 i-U-o.l i - 47 -0.05 i - 47 -0.05 

PL 2 norm 1.4±g;i x 10~ 4 5.0t{J;l x 10~ 4 S.6t° Q i x lO" 5 4.2±g:J x 10~ 4 2.4±g;| x 10~ 3 2.2± 1 ; 4 x 10~ 3 

Flux 1-10 keV l.lto! x lO" 10 l-l±g-| x 10- 10 2.6±gl x lO -11 2.6±J-| x 10- lx 2.2±g-| x 1Q- I1 2.1±gg x 10" 11 
Flux 1-200 keV2.3±\-l x 10- 10 2.3t 4 ;° x lO" 10 l.l±8;i x 10- 10 1.4±°j x 10- 10 l.l±g;| x lO" 10 1.1±8:| x lO" 10 

xl (dof) 0.99 (216) 0.80 (216) 1.11 (202) 1.01 (202) 1.14 (158) 1.08 (156) 



Note. — Best fit values of the spectral parameters obtained by fitting the ~l-200keV XMM-Newton and INTEGRAL AXPs' spectra with 
a blackbody plus two power-laws model (BB+2PL) for 4U Of 42+6f4 and f RXS J1708-4009 , while a single power-law was used for IE f 841-045 . 
Furthermore, all the sources were modeled with a resonant cyclotron scattering model plus a power-law (RCS+PL). Errors are at ftr confidence 
level, reported fluxes are absorbed and in units of ergs _1 cm — 2 , and in units of 10 22 cm -2 and assuming solar abundances from Lodders 
(2003); 2% systematic error has been included. See also Fig. [2] and S I4.1l for details. *: source slightly variable in flux and spectrum, see text 
for details. 
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4U 1 42+6 1 4 




Energy (keV) Energy (krrV) 



Fig. 2.— 4U 0142+614, 1RXS J 1708-4009 and IE 1841-045 : left column shows the spectra in Counts/s/keV while in the right column 
we report the v¥ v plots. For 4U 0142+614 and 1RXS J1708-4009 the upper panels are relative to the modeling with a blackbody plus two 
power-laws (BB+2PL), while we used a blackbody plus power-law for IE 1841-045 . Bottom panels report for all the sources the resonant 
cyclotron scattering plus a power-law model (RCS+PL). See Tab. [2] and § 14, II for details. 
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3 5 10 | 2 5 10 

Energy (keV) Energy (keV) 



Fig. 3. — IE 2259+586: left column shows the spectra in Counts/s/keV while in the right column we report the zaF„ plots. The upper 
panels are relative to the modeling with a blackbody plus two power-laws (BB+2PL), while bottom panels report the resonant cyclotron 
scattering plus a power-law model (RCS+PL). Note the hard X-ray spectrum has been fixed at the value from Kuiper et al. (2006). See 
Tab.l3land oETTlfor details. 



Table 3 

Spectral Parameters: IE 2259+586 



AXP 




IE 2259- 


-586* 




Parameters 


BB+ 


2PL 


RCS 


fPL 


N H 


0.97 


fO.04 

-0.03 


0.89^ 


-0.02 
-0.02 



kT (keV) 0.41+°;°! 0.32lg;{g 



BB norm 2.771^ x 10 



+0.02 x 1( ,-4 
0.03 



-0.02 



Ti 3.98 

PLi norm 4.89±^4 x'lO -2 

Pt 0.32t° ;U 3 

T 1 0+°- 2 

' res u — 0.2 

RCS norm l.O^ j x 10" 3 



r 2 1.02 1.02 

PL 2 norm 1.65±^ x 10~ 7 5.0±\^ x 10~ 5 

Flux 1-10 keV 2.5+^ x 10~ n 2.5lg;J x 10" 11 

xl (dof) 1.15 (178) 0.94 (178) 



Note. — Best fit values of the spectral parameters obtained by fitting the ~l-10keV XMM-Newton observation of fE 2259+586 with a 
blackbody plus two power-laws model (BB+2PL), and with a resonant cyclotron scattering model plus a power-law (RCS+PL). We fixed the 
second power-law photon index to Y2 = f .02, the value reported in Kuiper et al. (2006) from KXTE measurements. Errors are at lcr confidence 
level, reported fluxes are absorbed and in units of ergs - 'cm -2 , and N# in units of 10 22 cm -2 and assuming solar abundances from Lodders 
(2003); 2% systematic error has been included. See also Fig.[3]and 6 H,ll for details. *: source variable in flux and spectrum, see text for details. 
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12 5 10 l 2 5 io 

Energy (keV) Energy (keV) 



Fig. 4. — IE 1048-5937 : left column represents the spectra in Counts/s/keV while in the right column we report the v¥ v plots. The 
upper panels are relative to the modeling with a blackbody plus one power-law (BB+PL), while bottom panels report the resonant cyclotron 
scattering model (RCS). See Tab. [4] and § 14.21 for details. Black, blue, and light-green colors are relative to observations taken in 2007, 2005 
and 2003, respectively. The red lines represent the total model, while the dashed lines are the single components. 









Table 4 












Spectral Parameters: IE 1048-5937 






AXP 






IE 1048-5937 








2003 


2005 




2007 


Parameters 


BB+PL 


RCS 


BB+PL 


RCS 


BB+PL 


RCS 




L - b6 -0.03 


0-98™ 


J - DD -0.04 


73+ l) - U4 

u - ' °-0.04 


i 7-1+U.U4 
1 - 1 i -0.03 


0.82™ 


kT (keV) 
BB norm 

PLi norm 


O.GSjJg 

Longs x 10- 4 

q qi +0.02 

°- oi -0.04 

l.lOtlH x 10- 2 


0.39+™ 4 


64+ 03 

u - D4 -0.04 

0-7±o"i x 10 " 4 

O 1 O + 0.03 

J - io -0.04 

0.7±g"i x io- 2 


44+ - 05 

u -^-0.04 


7S+ - 01 
3.00+ u 8 x 10- 4 

o on+0 07 
J - zu -0.07 

2.2l°; 4 x IO" 2 


0.4518$ 


/3 T 

T~res 

RCS norm 




29+ 02 

u ' ZM -0.02 
o 7+0.2 
z -'-0.4 

1.9±° -\ x 10" 


u - oo -0.04 

2 Q+o.i 

- 4 i.oi±8:?? x io- 4 




29+ 05 
4 7+0-2 

' -0.2 

3.0l°; 4 x 10- 4 


Flux l-10keV l.lloi x 10 ~ n 


i.i±o;t x io- 


11 0.8±g;| x 10- 11 Oi 


^; 4 4 x io- 11 


3.0t° Q i x 10- 11 


3.o±8:2 x io- 11 


xl (dof) 


0.99 (176) 


0.98 (176) 


0.99 (153) 


1.00 (153) 


1.08 (184) 


1.23 (184) 



Note. — Best fit values of the spectral parameters obtained by fitting several ~l-10keV XMM-Newton spectra, taken in different source 
states, with a blackbody plus power-law model (BB+PL), and with a resonant cyclotron scattering model (RCS). Errors are at la confidence 
level, reported fluxes are absorbed and in units of ergs - 1 cm — 2 , and Nj^ in units of 10 22 cm -2 and assuming solar abundances from Lodders 
(2003); 2% systematic error has been included. See also Fig. [4] and § 14.21 for details. 
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XTE 1B10-197 



1; 




XTE 1810-197 




XTE 1B10-197 




XTE 1B10-197 




Energy (keV) Energy (keV) 

Fig. 5. — XTE J1810-197 : left column represents the spectra in Counts/s/keV while in the right column we report the uF v plots. The upper 
panels are relative to the modeling with two absorbed blackbodies (BB+BB), while bottom panels report the resonant cyclotron scattering 
model (RCS). See also Tab. [3] and § 14,21 for details. Black, light-green and blue colors are relative to observations taken on 2004, 2005 and 
2006, respectively. The red lines represent the total model, while the dashed lines are the single components. 



Table 5 

Spectral Parameters: XTEJ1810-197 



AXP 



XTEJ1810-197 



2004 



2005 



2006 



Parameters 



BB- 



BB 



RCS 

TTJW 



BB- 



BB 

-0.07 



RCS 

OTDTT 



BB+BB 



RCS 
for - 



N 



H 



0.58; 



0.40 



-0.05 



0.52 



0.25 



-0.05 



0.14 



kTi (keV) 
BBi norm 
kT 2 (keV) 
BB2 norm 



6.6 



0.36] 
+0.5 
-0.4 

0.71 



12 + n x 10" 



-0.02 
x 10" 
+0.01 
-0.02 



0.44 



-0.03 
-0.03 



0.27: 



hO.03 
-0.02 

3.8±g-? x 10~ 5 

u -°°-0.03 

i.5tn x 10- 5 



0.291°°* 



0.25; 



2 7 +0 - 3 x 

z -'-0.3 x 



0.04 

10" 



^fi+ 05 

u - do -0.07 



0.7 



+0.1 
-0.1 



x 10" 



0.13 



-0.05 
-0.05 



0T 
7~res 

RCS norm 



1.1" 



+0.05 
0.05 
+ 1.6 
-1.0 

0.3 >< 10- 



0.19 

5.9 



40+ 05 

u - w -0.05 

1 6 +0 - 2 
7.2t° i x 10~ 5 



S5+ 05 

u - oo -0.05 



1.4 

9 £+0.5 
z -°-0.5 



+0.1 
-0. 

X 



10" 



Flux 1-10 keV 12±l 



x 10" 



x 10" 



2.2 



-0.1 



x 10" 12 2.1 



hO.2 



x 10" 12 1.2 



+0.3 



x 10" 12 1.2 



-0.5 



x 10" 



Xl (dof) 



1.21 (135) 1.27 (135) 0.94 (97) 



1.07 (97) 



0.97 (67) 



1.00 (67) 



Note. — Best fit values of the spectral parameters obtained by fitting several ~1— 10 keV XMM-Newton spectra, taken in different source 
states, with two absorbed blackbodies (BB+BB), and with a resonant cyclotron scattering model (RCS). Errors are at Icr confidence level, 
reported fluxes are absorbed and in units of ergs - 'cm" 2 , and Nj^ in units of 10 22 cm -2 and assuming solar abundances from Lodders (2003); 
2% systematic error has been included. See also Fig. [5] and § 14.21 for details. 
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12 5 10 i 2 5 10 

Energy (kcVl Energy (keV) 



Fig. 6. — IE 1547.0-5408 : left column represents the spectra in Counts/s/keV while in the right column we report the vF v plots. The 
upper panels are relative to the modeling with two blackbodies (BB+BB), while bottom panels report the resonant cyclotron scattering model 
(RCS). See also Tab. [5] and § 14,21 for details. Black and light-green colors are relative to observations taken on 2007 and 2006, respectively. 
The red lines represent the total model, while the dashed lines are the single components. 



Table 6 

Spectral Parameters: IE 1547.0-5408 



AXP IE 1547.0-5408 

2006 2007 

Parameters BB+BB RCS BB+BB RCS 

3J6™ 2^ 458™ 46l£T 



k Tl (keV) 0.46±g:g 0.331°;°* 0.51±g : g 0.46_t°;°| 

BBi norm 1.2±g"| x 10~ 5 7.2±g-| x 10~ 6 

kT 2 (keV) 1.2±g:i l-34ig;g? 

BB 2 norm 1.4±°;^ x 10~ 6 l-4±o.i x 10~ 4 



a n QO+0.03 n O/ ,+0.04 

Pt U.32_ 09 U.24_ 004 

1 o +0 - 8 1 0+ ' 1 

'res ± - u — 0.2 ± ' u — 1 

RCS norm 2.6±g;| x 10~ 5 9.4_t° j x 10" 



Flux l-10keV3.2±°^ x 10- 13 3.1+°; 2 x 10" 13 3.0l^ x 10~ 12 3.01°^ x 10~ 12 
xl (dof) 1.18 (60) 1.20 (60) 1.02 (105) 1.13 (105) 



Note. — Best fit values of the spectral parameters obtained by fitting several ~1— 10 keV XMM-Newton spectra, taken in different source 
states, with two absorbed blackbodies (BB+BB), and with a resonant cyclotron scattering model (RCS). Errors are at Icr confidence level, 
reported fluxes are absorbed and in units of ergs - 'cm" 2 , and Nj^ in units of 10 22 cm -2 and assuming solar abundances from Lodders (2003); 
2% systematic error has been included. See also Fig. [6] and § 14.21 for details. 
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Fig. 7. — CXOU J1647-4552 : left column represents the spectra in Counts/s/keV while in the right column we report the i/F v plots. 
The upper panels are relative to the modeling with two absorbed blackbodies (BB+BB), while bottom panels report the resonant cyclotron 
scattering model (RCS). See also Tab.[7]and £j |4.2l for details. Black and light-green colors are relative to observations taken on 2006 September 
22 and 16, respectively. The red lines represent the total model, while the dashed lines are the single components. 



Table 7 

Spectral Parameters: CXOU J1647-4552 



AXP 



CXOUJ1647-4552 



2006/09/16" 



2006/09/22" 



Parameters 



BB- 



BB 



RCS 

TTTET 



BB- 



BB 



RCS 



N 



H 



2.14 



2.08 



-0.16 



2.34" 



2.40 



kTi (keV) 
BBi norm 
kT 2 (keV) 
BB2 norm 



0.39 



+0.03 
0.02 



0.85± i 

2 - 4 -8:i x 10-6 



0.34; 



-0.15 
-0.19 



0.59: 

1-0.5 



0.02 



1.4 



1.23 

+0.1 



+0.04 
-0.04 

x 10" 



0.55 



+0.08 
-0.08 



Pt 

7~res 

RCS norm 



0.30 
2.9 



0.08 
+0.1 



7.8to.l x 10" 



0.08 
0.08 
+0.05 
0.05 

3.0t°oi x 10" 



0.42 
1.09 



Flux 1-10 keV 2-4±g;i x 10- 13 2.4+^ x 10" 13 2.2±°;^ x 10~ n 2.2±%\ x 10~ n 
xl (dof) 1.00 (73) 1.23 (73) 1.01 (136) 1.06 (136) 



Note. — Best fit values of the spectral parameters obtained by fitting several ~1— 10 keV XMM-Newton spectra, taken in different source 
states, with two absorbed blackbodies (BB+BB), and with a resonant cyclotron scattering model (RCS). Errors are at Icr confidence level, 
reported fluxes are absorbed and in units of ergs - 'cm" 2 , and Nj^ in units of 10 22 cm -2 and assuming solar abundances from Lodders (2003); 
2% systematic error has been included. See also Fig. [7] and § 14.21 for details. 
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Fig. 8. — SGR 1806-20 : left column shows the spectra in Counts/s/keV while in the right column we report the jaF„ plots. The upper 
panels are relative to the modeling with a blackbody plus power-law (BB+PL), while bottom panels report the resonant cyclotron scattering 
model plus power-law (RCS+PL). See also Tab. [8] and £14.31 for details. Light green, black and blue colours are relative to observations taken 
on 2003, 2004 and 2005, respectively. The red lines represent the total model, while the dashed lines are the single components. 



Table 8 

Spectral Parameters: SGR 1806-20 



SGR SGR 1806-20 



2003 2004 2005 



Parameters BB+PL RCS+PL BB+PL RCS+PL BB+PL RCS+PL 

~n~h 9^ 9^ WTY: ul2-^ ii.o,:;! 

kT(keV) 0.56™ 0.57±°;f 0.72™ 0.54™ .57t°;°* OJ86±gS 

BB norm 5.5±^J x 10~ 5 l-0±a3 x 10 ~ 4 1A -H x 10 ~ 5 

a n -i 7 +0.03 n lzL +0.08 n ZLQ+ - 04 

PT U.l/_ 03 U.14_ 03 0.49_ 003 

2 2+ 1 ' 5 4 3 +0 - 7 2 6 +0 - 2 

I res i ^ t -'->_l i ^"U—O 3 

RCS norm 3.8±g;| xlO" 5 7.4^] x lO" 5 4.6±° l x'l0~ 4 

r 1 x+0.1 i o+O.l -, o+0.1 -i o+O.l -, r+0.1 i o +0.2 

1 i -°-o.i i - z -o.i 1 -°-o.i i - J -o.i J -°-0.1 i - z -o.i 

PL norm 3.1±g| x 10~ 3 l-7±g;| x 10~ 3 4.7±gj x 10 ~ 3 S.lt^j x 10~ 3 3.8±g;l x 10" 3 1.7±^ x 10~ 3 

Flux l-10keV1.2±{J;jj x 10 -11 1.2±g;l x 1Q- 11 2.6t° 6 7 x lO" 11 2.6±g:| x 10- 11 1.4±^;| x lO" 11 l-3±g;l x lO" 11 

xl (dof) 0.96 (54) 1.03 (52) 1.01 (65) 0.97 (63) 1.02 (159) 0.90 (157) 



Note. — Best fit values of the spectral parameters obtained by fitting several ~l-10keV XMM-Newton spectra, taken in different source 
states, with a blackbody plus power-law model (BB+PL), and with a resonant cyclotron scattering plus power-law model (RCS+PL). Errors 
are at lcr confidence level, reported fluxes are absorbed and in units of ergs -1 cm — 2 , and Njj in units of 10 22 cm - 2 and assuming solar 
abundances from Lodders (2003); 2% systematic error has been included. See also Fig. [8] and i|4.3l for details. 
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Fig. 9. — SGR 1900+14: left column shows the spectra in Counts/s/keV while in the right column we report the z/F„ plots. The upper 
panels are relative to the modeling with a blackbody pl us p ower-law (BB+PL), while bottom panels report the resonant cyclotron scattering 
model plus power-law (RCS+PL). See also Tab. [9] and £|4.3I for details. The red lines represent the total model, while the dashed lines are the 
single components. 



Table 9 

Spectral Parameters: SGR 1900+14 



SGR 


SGR 1900+14 


Parameters 


BB+PL RCS+PL 


N H 


o c+U.i A O+U.l 

O.O_ x 4 - u -0.1 


constant 


1.20 1.10 



kT (keV) 0.45±^ 0.30^'f 

o.i 
o.i 



BB norm 6.712 i x 10 



(b 0.26t°;°£ 

t 2 5 +0 ' 5 

' res ^•" J _o 2 

RCS norm l-8lo.05 x 10 ~ 4 

T 1 4 +al 1 24+ 07 

PL norm 4.4l°;J x 10" 4 3.01°;* x 10~ 4 



Flux 1-10 keV 3.9ljH x 10- 12 3.8l^ x 10 
Flux 1-200 keV1.7lg;i x 10" 11 1.7±o\ x 10 



-12 
11 



xl (dof) 1.18 (141) 1.15 (139) 



Note. — Best fit values of the spectral parameters obtained by fitting the ~l-200keV XMM-Newton and INTEGRAL spectra with a 
blackbody plus a power-law model (BB+PL), and with a resonant cyclotron scattering model plus a power-law (RCS+PL). Errors are at la 
confidence level, reported fluxes are absorbed and in units of crgs _1 cm~ 2 , and Nj/ in units of 10 22 cm -2 and assuming solar abundances from 
Lodders (2003); 2% systematic error has been included. See also Fig. [9] and i|4,3l for details. 
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Fig. 10. — Angle- averaged optical depth in a twisted magnctosphere model (Thompson Lyutikov & Kulkarni 2002) as a function of the 
twist angle. The curve refers to Ptulk = 1; f° r different values of the bulk velocity the ordinate should be divided by fibulk- 
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Fig. 11. — Comparison between the derived spectral parameters and the sources' properties (see [[5]for details). To infer the l-10keV 
luminosity we assumed a distance of 3, 3, 5, 3.3, 3, 7, 10, 4, 5, and lOkpc, for the sources ordered as the labels reported in the top panel 
(from left to right and top to bottom). Errors in the luminosities are assumed to be 30% of the reported values (which is of the order of the 
flux errors), although the real error (including that on the distance) is actually much larger. 



